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I.  INTRODUCTION 

This  report  summarizes  the  work  performed  during  the  contract 
period,  July  1,  1974  through  November  30,  1977.  It  should  be  noted 
a number  of  the  studies  are  being  continued  under  a new  ARO  contract 
which  commenced  on  January  1,  1978. 


II.  DESCRIPTION  OF  PROBLEMS  STUDIED 


The  basic  objective  of  the  project  was  to  develop  and  apply  the 
technique  of  computer  molecular  dynamics  simulation  to  problems  in 
materials  science  research.  Computer  programs  were  developed  for 
the  simulation  of  atomic  and  molecular  solids  in  which  atoms  interact 
through  continuous  interatomic  potentials.  Computer  programs  also 
have  been  developed  for  the  simulation  of  hard-sphere  mixtures  which 
undergo  chemical  reactions.  Applications  of  these  programs  have 
been  made  to  study  the  structure  and  motions  of  grain-boundary 
systems,  molecular  vibrations  and  reorientations  in  solids,  and  in- 
stabilities in  multicomponent  systems  undergoing  reactions. 

* 

Equilibrium  Properties  and  Dynamics  of  Grain  Boundaries 

A basic  problem  in  the  study  of  thermodynamic  properties  of 
grain-boundary  systems  is  the  calculation  of  the  boundary  surface 
tension. ^ At  a given  temperature  the  surface  tension  can  be  ex- 
pressed in  terms  of  the  excess  internal  energy  and  the  excess  entropy 
(relative  to  the  reference  system  which  is  the  perfect  solid) . We 
have  developed  a method  for  computing  the  vibrational  part  of  the 
entropy  by  using  computer  simulation  to  generate  the  vibrational 

frequency  spectrum  of  the  boundary  system  and  to  evaluate  the  system 

(2) 

internal  energy  over  a range  of  temperature.  The  frequency 
spectrum  is  used  to  calculate  the  entropy  in  the  quasiharmonic  ap- 
proximation at  a sufficiently  low  temperature  where  the  approximation 
is  valid,  and  the  additional  entropy  is  obtained  by  integrating  the 
ratio  of  the  internal  energy  to  the  temperature  from  this  low  tem- 
perature up  to  the  temperature  of  interest.  This  method  has  been 
applied  to  a 400-atom  three-dimensional  system  containing  two  large- 
angle  tilt  grain  boundaries  and  a structural  vacancy. 


This  part  of  the  project  has  been  carried  out  in  collaboration  with 
Drs.  G.  Bishop,  A.  Cox,  and  R.  Harrison  of  Army  Materials  and 
Mechanics  Research  Center,  Watertown,  Mass. 
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Motions  of  grain  boundaries  are  of  fundamental  importance  in 
understanding  the  behavior  of  materials  at  elevated  temperatures. 

The  coupled  motion  of  boundary  sliding  and  migration  has  been  dis- 
cussed in  the  literature  on  the  basis  of  theoretical  considerations 

(3) 

and  bubble-raft  models.  We  have  observed  this  process  in  the 

(4) 

course  of  our  simulation  of  grain-boundary  systems.  A series  of 
simulation  experiments  were  performed  to  study  the  effects  of  tem- 
perature, initial  stress  on  the  system,  number  of  atoms  in  the  system, 
and  the  type  of  interatomic  potential  function  used.  In  addition, 
simulations  were  carried  out  on  two-dimensional  grain-boundary  systems 
which  do  not  contain  any  structural  vacancies  and  recrystallization 
of  the  solid  was  observed  as  a competing  process  to  boundary  sliding 
and  migration. 

Molecular  Solids 

We  have  developed  a simulation  program  which  treats  a molecular 

solid  as  a system  of  atoms  which  interact  through  bonding  and  non- 

bonding  potential  functions. This  program  evolved  from  an  earlier 

version  of  simulation  program  for  atomic  solids  which  was  used  to 

compute  the  thermodynamic  and  lattice  dynamic  properties  of  solid 

argon  over  a wide  range  of  temperatures  and  pressures. ^ Simulation 

of  molecular  solids  thus  far  has  been  confined  to  a-^  where  time- 

dependent  correlation  function^  describing  molecular  vibrations  and 

reorientations  have  been  computed. ^ Recently  a study  was  initiated 

to  investigate  the  vibrational  relaxation  of  a diatomic  molecule  in 

an  atomic  lattice. ^ it  is  expected  that  this  work  will  lead  to  a 

study  of  dissociation/recombination  dynamics,  in  particular,  the 

(8) 

problem  of  geminate  recombination  in  gases  and  liquids. 

Chemical  Reacting  Systems 

The  technique  of  computer  molecular  dynamics  largely  has  not 
been  applied  to  problems  of  chemical  reactions  thus  far.  We  have 
performed  a number  of  exploratory  calculations  on  hard-sphere  systems 
which  show  that  the  important  features  such  as  interspecies  conversion. 


IT 
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spontaneous  decay,  and  release  of  internal  energy  of  a molecule  during 
collisions  all  can  be  simulated  in  a computer  molecular  dynamics  ex- 
periment. In  particular  we  have  obtained  computer  results  useful  for 

(9) 

testing  the  validity  of  conventional  macroscopic  rate  equations. 

The  computer  program  used  in  this  work  is  an  adaptation  of  a program^1 
written  originally  to  study  hydrodynamic  flow  processes  and  was  ef- 
fectively a program  for  two-dimensional  systems.  We  have  recently 
generated  a corresponding  program  for  three-dimensional  systems  with 
periodic  boundary  conditions.  One  of  the  applications  of  this  program 
will  be  the  study  of  collision- induced  de-excitation  of  internal 
energy  states  as  a model  of  the  process  of  thermal  explosion.  Other 
problems  to  be  investigated  in  the  continuation  of  this  work  include 
chemical  Instabilities  in  multicomponent  systems  and  dynamical  proper- 
ties of  diffusion-reaction  systems. 
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III.  SUMMARY  OF  MOST  IMPORTANT  RESULTS 

The  most  significant  result  to  emerge  from  our  research  are  the 
observations  of  grain  boundary  sliding  and  migration  in  finite- 
temperature  simulation  of  three-  and  two-dimensional  solids.  These 
results  are  the  first  demonstration  that  grain-boundary  motions  can  be 
simulated  under  fully  dynamic  conditions.  They  also  serve  to  confirm 
the  validity  of  the  concept  of  "DSC  lattice"  vectors  of  coincidence 
site  lattice  theory.  Detailed  analysis  of  the  simulation  results 
provide  considerable  insight  Into  the  effects  of  the  local  strain 
generated  by  the  defect  in  determining  the  behavior  of  grain  boundaries. 

Another  significant  result  of  our  study  of  grain-boundary  systems 
is  the  method  for  computing  vibrational  entropy  through  the  vibrational 
frequency  spectrum  without  the  assumption  of  harmonic  forces.  By 
using  computer  simulation  data  in  conjunction  with  the  techniques  of 
lattice  dynamics  one  can  considerably  extend  our  capability  of  analyz- 
ing defect  structures  in  solids  at  elevated  temperatures. 

In  the  area  of  chemical  reactions  we  have  provided  the  first 
demonstration  that  reactions  involving  interspecies  conversion  and 
spontaneous  decay  can  be  studied  In  detail  by  computer  molecular 
dynamics  simulation.  We  have  obtained  results  which  showed  that 
multiple  steady-state  solutions  predicted  by  the  analysis  of  macro- 
scopic rate  equations  indeed  occur  and  that  the  rate  equations  give  a 
good  description  of  the  average  behavior  of  the  reacting  system.  We 
have  also  demonstrated  the  feasibility  of  simulating  exothermic 
reactions  with  thresholds,  thus  providing  the  basis  for  studying  energy 
dissipation  mechanisms  in  highly  unstable  systems. 
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